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1 Introduction 



The largest part of the cross section for inelastic processes at LEP2 is due to two-photon 
scattering. Interactions of two photons can be studied at e + e~ colliders by investigating 
the reaction 

e + e~ — > e + e~77 — > e + e~X (1) 

where the photons can be quasi-real or virtual and the hadronic final state is denoted by 
X. 

The analysis reported in this paper focusses on the interactions of virtual photons 
by selecting only double tagged events, i.e., events where both scattered electrons are 
detected. Here, electrons and positrons are generically referred to as electrons. Differential 
cross sections for this process are measured as a function of several kinematic variables 
and presented here. 

The purpose of this paper is to compare the measured differential cross sections with 
the predictions of the PYTHIA [L and PHOT02 [2] Monte Carlo generators, with a next- 
to-leading-order (NLO) QCD calculation [3] and with BFKL calculations HE1E10IB! • 

2 Theoretical Framework 

Interactions of virtual photons can be studied by requiring that both electrons be detected 
after radiating the photons. The kinematics of these electron-induced 77 interactions 
is sketched in Fig. ^ The symbols in parentheses represent the four-momenta of the 
particles. 




Figure 1: Kinematics of 77 interactions at an e + e collider 
The kinematics can be described by the dimensionless Bjorken variables of deep 
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inelastic scattering: 



T . = 91 (0) 

Qj + Ql + W^ ' [Z) 

Vi = I-t^cos 2 ^, (3) 



Eh 



beam 



where i — 1,2 refers to the scattered electron, positron. The hadronic invariant mass W 77 
used in these definitions is obtained from the energies Eh and the momenta ph of the final 
state particles h by 

K = (e ^) - (l>) = E l - • ( 4 ) 

The virtualities of the photons Q\ are 

Qi = -iPi - P'if = 2E hcam E[(l - cos 6i) (5) 

for 6i 3> m e i ec t r on/ ' Ei, where 6i are the scattering angles of the deflected leptons. 

From Eq. (j3J) it is possible to select interactions of virtual photons by requiring that 
the scattered electrons be detected at large angles. The accessible range of virtualities 
and therefore the phase space depends on the region where electrons can be detected. 

For the comparison of the data with BFKL calculations the following quantity is 
defined, 

Y = log (-) , (6) 



where y/s is the LEP centre-of-mass energy and s/sq = syxy-ij \JQ\Q\ ~ ^ 77 / \jQ\Q\i 



the approximation requiring Wjr 3> Q~. 

To correct the data for detector effects, the signal Monte Carlo simulation for this 
analysis uses PYTHIA 6.151 [UHEEm and as an alternative PHOT02 0. The PHOJET 
generator [11] . which simulates untagged events, gives information for further systematic 
studies. 

The PHOT02 generator is a combination of different two-photon generators. The 
main contribution to this analysis comes from the QED part, based on a program 
that gives a matrix element calculation for e + e~ — > e + e~£ + £~, where £ ± are charged 
fermions [T2* | ll3 [ ll4]. A small contribution arises from the vector- meson dominance model 
( VDM) [THjlTT)] . The implementation of the 77 generator in PYTHIA is based on a model 
described in pHITIlj. 

The background from leptonic double tagged two-photon events is also simulated 
with PHOT02. For annihilation events three generators are employed: KORALZ [T7] . 
PYTHIA and HERWIG jTHJ. The latter is used to generate single tagged two-photon 
events contributing to the background. 
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3 ALEPH Detector 



A detailed description of the ALEPH detector and its performance can be found in 
Ref. The inner part of the ALEPH detector is dedicated to the reconstruction 

of the trajectories of charged particles with a two-layer silicon strip vertex detector 
(VDET), a cylindrical drift chamber (ITC) and a large time projection chamber (TPC). 
The three tracking detectors are immersed in a 1.5 T axial magnetic field provided by a 
superconducting solenoidal coil. Together they measure charged particle momenta with 
a resolution of Sp t /pt = 6 x 10~ 4 p t © 0.005 (p t in GeV/c). 

Photons are identified in the electromagnetic calorimeter (ECAL), situated between 
the TPC and the coil. The ECAL is a lead/proportional-tube sampling calorimeter 
segmented into 0.9° x 0.9° projective towers read out in three sections in depth. It 
has a total thickness of 22 radiation lengths and yields a relative energy resolution of 
0.18/ \/E + 0.009 (E in GeV) for isolated photons. Electrons crossing the TPC are 
identified by their transverse and longitudinal shower profiles in ECAL and their specific 
ionization in the TPC. 

The iron return yoke is instrumented with 23 layers of streamer tubes and forms the 
hadron calorimeter (HCAL). The latter provides a relative energy resolution for charged 
and neutral hadrons of 0.85/ \[E [E in GeV). Muons are distinguished from hadrons by 
their characteristic pattern in HCAL and by the muon chambers, which are composed of 
two double-layers of streamer tubes outside HCAL. 

The information from the tracking detectors and the calorimeters are combined in an 
energy-flow algorithm [T^]. For each event, the algorithm provides a set of charged and 
neutral reconstructed particles, called energy-flow objects. 

Two small-angle luminosity calorimeters, the silicon luminosity calorimeter (SICAL) 
and the luminosity calorimeter (LCAL), are used to detect and measure the energies 
of the electrons from beam-beam scattering including the electrons in the final state of 
reaction (1). The SICAL uses 12 silicon/tungsten layers to sample showers. It is mounted 
around the beam pipe and covers angles from 34 to 58mrad. An energy resolution of 
0.225-^^ (E in GeV) is achieved. The LCAL is a lead/proportional-tube calorimeter, 
similar to ECAL, placed around the beam pipe at each end of the ALEPH detector. It 
monitors angles from 45 to 160mrad with an energy resolution of 0.33v^E (E in GeV). 

4 Data Sample 
4.1 Event Selection 

The analysis is based on the data taken with the ALEPH detector from 1998 to 2000 at 
centre-of-mass energies \fs = 188 — 209 GeV and corresponds to an integrated luminosity 
of 640 pb -1 . 

The event selection is performed in three stages: detection of scattered electrons, 
verification of the presence of a hadronic system, and background reduction. 
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• Detection of scattered electrons. 

The luminosity detectors SICAL and LCAL detect the scattered electrons. Thus 
the polar angular range is restricted to 35 mrad < 9i < 55 mrad (SICAL) and 
60mrad < Q{ < 155 mrad (LCAL). The energy threshold is set to E[ > 0.3.E'beam- 

• Verification of the hadronic system. 

To ensure that the final state is a hadronic system and not a lepton pair, at least 
three charged particles are required. The visible mass W 77 of the hadronic system 
must be larger than 3GeV/c 2 . 

• Background reduction. 

The total visible energy E to t = E[ + E' 2 + Ex must be larger than 70% of the nominal 
centre-of-mass energy. To reject remaining Bhabha events the acolinearity of the 
scattered leptons is required to be less than 179.5°. 

4.2 Backgrounds 

With these cuts 891 events were selected in the data with 206.1 expected background 
events. The three remaining sources of background are the following. 

• Double-tagged leptonic 77 events containing mostly r + r _ as estimated with the 
PHOT02 Monte Carlo simulation. 

• Superpositions of single tagged 77 events and off-momentum electrons. In order to 
appraise this background source it is necessary to extract the probability of finding 
an off- momentum electron in an arbitrary data event. This is done by looking for 
additional energy deposits from off-momentum electrons in Bhabha events. These 
off-momentum electrons are then added to single tagged events simulated with the 
HERWIG Monte Carlo or alternatively to those taken from data. 

• Annihilation events (e + e~ — > qq) which can also fake the topology of double tagged 
events. This background source, dominated by radiative Z return, is generated using 
the KORALZ Monte Carlo. 

The trigger efficiency is estimated by comparing the rates of two independent triggers: 
the Bhabha event trigger and the combination of non-Bhabha (charged-track and neutral- 
energy) triggers. The selected events are found to allways fulfill the two trigger conditions. 
Therefore, the trigger is taken to be 100% with an uncertainty of 2% and no correction is 
applied. 

4.3 Selection Results 

The numbers of events obtained in data, the signal Monte Carlo simulations and the 
estimated backgrounds are summarized in Table H The total cross section of the signal 
Monte Carlo simulations was normalized to data after background subtraction. The cross 
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Table 1: Numbers of observed events for combinations of detectors. The notation "MC + back" stands 
for PYTHIA Monte Carlo plus all background sources. PHOT02 is given for comparison as an alternative 
signal Monte Carlo. 



SiCAL - SiCAL 


243 


195 


150 


177 


13 


3 


30 


SiCAL - LCAL 


388 


447 


360 


328 


37 


30 


21 


LCAL - LCAL 


260 


249 


176 


179 


23 


51 





total 


891 


891 


685 


685 


73 


83 


50 



section of the PYTHIA Monte Carlo generator was reduced by 12% and the cross section 
extracted from PHOT02 was increased by 30%. 

Several measured spectra are given in Figs. 121 and El showing the comparison between 
data and simulations after normalization. 



5 Acceptance Corrections and Systematic Errors 

A simple bin-by-bin method was applied to correct for detector inefficiencies. The 
correction factors were calculated for each bin as 

Atrue ,„\ 

Shin = -7Z ( 7 ) 

'detected 

where Atme is the number of generated events in a given bin and Adetected the number of 
detected events in the same bin according to the simulation. The corrected data values 
for a bin, -R cor , were then derived from the measured values, invisible, by 



-*^cor (,' '-visible -^background ) ^bin (8 J 

where -Rbackground is the expected background. 

In order to estimate the systematic effect caused by an imperfect detector simulation 
all energy resolutions were varied by 10% (Fig. EJ). 

The uncertainty due to a possible shift in the energy scale of the luminosity monitors 
SICAL and LCAL (Fig. HJ was estimated by introducing an offset of 0.5 GeV to the 
measured energy. The polar and azimuthal angles of the electrons were shifted by 
0.25 mrad and 0.5mrad. These seemingly large adjustments were made to estimate a 
possible systematic uncertainty due to a poor description of small polar angles of the 
electrons (Fig. [H}. 

The cross sections of the background processes were changed conservatively by ±10%. 
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The PHOT02 Monte Carlo simulation was used instead of PYTHIA to correct for 
detector effects. In both cases the systematic uncertainties due to statistical fluctuations 
in the Monte Carlo sample were small. 

Systematic differences between the data collected in different years were not observed. 

The systematic error was computed as the quadratic sum of the various contributions. 
The main contributions to the systematic error come from varying the energy resolutions 
and the tag energy. 

6 Results 

The measured cross sections, corrected for detector effects, are given in Figs. E]to[T3] The 
corresponding bin-contents are given in the Appendix. The results are compared with the 
PYTHIA and PHOT02 Monte Carlo models and the NLO QCD prediction 0. 

Figure |U shows the cross section as a function of the tag energy. The cross section 
predicted by the NLO QCD calculation comes out too low by 20%. The tail towards 
low energies is slightly underestimated by all models. The polar angle 9 of the scattered 
electrons is given in Fig. The first measured point (35 to 40 mrad) is overestimated by 
the Monte Carlo. The gap between SICAL and LCAL (55 to 60 mrad) is interpolated for 
the further results in Figs. EH to fl3l 

The virtuality of the photons Q\ are well simulated over the whole accessible range 
(Fig.EJ). This also applies to the ratio of the two virtualities (Fig. EJ). 

The acoplanarity angle A0 (Fig. |HJ) and the acolinearity $ between the scattered 
electrons (Fig. |HJ) are well described by the PYTHIA prediction. The NLO QCD 
calculation yields a slightly too low cross section at large angles, and PHOT02 does 
not describe well the acoplanarity angle distribution. The mass of the hadronic system is 
again well reproduced by the Monte Carlo simulations. The NLO QCD calculation comes 
out with a slightly too low cross section for small masses (Fig. [TO]) . The deep inelastic 
scattering variables Xi (Fig. ITT]) and yt (Fig. IT2"]) in data are well described by the Monte 
Carlo simulations. The NLO QCD calculation fails in reproducing the cross section as a 
function of Xi, while yi agress well. The cross section as a function of Y is given in Fig.fTTfl 
Again the Monte Carlo simulations describe the measured spectrum well. 

Finally the data are compared to BFKL predictions. The calculation was done for 
photons with virtualities of 12 GeV 2 (90% of the data) and 38 GeV 2 (10% of the data). 
Since this calculation assumes that both photons have the same virtuality, a further cut 
was added to the event selection [2TJj : 

\AQ\ = \\ogQ 2 1 /Ql\<1.0 . (9) 

The whole analysis, including systematic error estimation, was redone with this new cut, 
which reduced the statistics by 40%. 

The <t 77 cross section j2JJ|221l2H] was ex t rac t e d from the measured e + e~ cross section 
using G ALUGA [21] . The resulting cross section as a function of Y is plotted in Fig. ITU 
where it is compared to LO BFKL and NLO BFKL calculations HElinilZllH] ■ In the LO 
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BFKL calculation the Regge scale parameter was varied from Q 2 to 10Q 2 . Even at the 
lowest Y values this calculation is barely in agreement with the data. The NLO BFKL 
calculation, however, with the Regge scale parameter varied from Q 2 to 4Q 2 (solid curves) 
is consistent with the data. 

7 Conclusions 

The cross section for the process 

e + e~ — > e + e~77 — > e + e~hadrons (10) 

has been measured using ALEPH data taken at e + e~ centre-of-mass energies yfs = 
188 — 209 GeV and corresponding to an integrated luminosity of 640pb _1 . The phase 
space is defined by the electron energies E[ 2 > 0.3-Ebeam, the polar angles of the electrons 
0.35 mrad < < 155mrad, and the mass of the hadronic system W 77 > 3GeV/c 2 . 

The differential cross section for e + e~ — > e + e~77 — > e + e~X, X being the hadronic final 
state, has been measured as a function of various event observables. 

The majority of measured distributions are well described in shape by the PHOT02 and 
PYTHIA Monte Carlo models, but both required an adjustment to their normalization 
to match the data. The cross section of the PYTHIA Monte Carlo generator was reduced 
by 12% and the cross section extracted from PHOT02 was increased by 30%. 

The NLO QCD prediction yields a slightly low cross section. With the exception of 
the deep inelastic scattering variable Xi this calculation also gives a reasonable description 
of the data. Similar conclusions have been drawn by the L3 j2Sj and the OPAL j2Ej 
Collaborations. 

A slight enhancement of the data with respect to the simulations is observed at high 
Y, but the steep rise of the cross section as predicted by LO BFKL is not observed in 
data. However, NLO BFKL is close in shape and normalization to the data. 
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ALEPH DATA 

— PYTHIA 

- - PHOT02 

yy ^ tt 
nil momentum 

qq 



-m- ALEPH DATA 
■-- PYTHIA 
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yy — > tt 
oft momentum 

qq 




0.1 0.125 0.15 

e,, 2 [deg] 



-m- ALEPH DATA 
--- PYTHIA 
- - PHOT02 

off momentum 

qq 




150 200 
Qj [GeV 2 ] 
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E2 yv^tt 
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AO [rad] 





-m- ALEPH DATA 




■■- PYTHIA 
- - PHOT02 
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[53 off momentum 
ID] qq 







1 




Figure 2: Differential cross sections as a function of various observables. Shown are (from top left 
to bottom right): The relative tag energy E[ 2 / '-Ebcam; the polar angle of the scattered electrons 01,2] 
the virtualities Q\ of the photons; the ratio of the two virtualities, shown as AQ = — logQi/Q^', the 
acoplanarity angle A<fi between the scattered electrons; the acolinearity $ between the scattered electrons; 
the deep inelastic scattering variables Xi and yi. The plots contain the background contributions (filled 
areas), the PYTHIA and PHOT02 Monte Carlo predictions plus background (lines), and the data (points) 
with the statistical errors. 
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-m- ALEPH DATA 

— PYTfflA 

- - PHOT02 

yy^TT 
oft' momentum 

qq 




-m- ALEPH DATA 
■-- PYTHIA 
- - PHOT02 

ES3 yy — > xr 

off momentum 
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Y=ln(W^/(Qfc2)" 2 ) 
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qq 
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qq 
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off momentum 

qq 
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S 0-4 



4 6 
^neutral 
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- - PHOT02 

\22 yy^Tx 

off momentum 
mi qq 



tot beam 



P, [GeV] 



Figure 3: Differential cross sections as a function of various observablcs. Shown are (from top left to 
bottom right): The mass of the hadronic system W 77 ; the quantity Y = log W^/y/ Q\Q\\ the number 
of charged tracks N c h\ the pseudorapidity r] c h = — logtan(f9 c h/2) of the charged tracks; the number of 
neutral objects N neutla \; the pseudorapidity 77 noutra i of the neutral objects; the relative total visible energy 
i?tot/2i?bcam; the total transverse momentum p t of the event with respect to the beam pipe. The plots 
contain the background contributions (filled areas), the PYTHIA and PHOT02 Monte Carlo predictions 
plus background (lines), and the data (points) with the statistical errors. 
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Figure 4: Differential cross section as a function of the energy of the scattered electrons Ei^/Eb cam 
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Figure 5: Differential cross section as a function of the polar angle 9 of the scattered electrons 
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Figure 6: Differential cross section as a function of the virtualities Q\ 2 of the two photons 
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Figure 7: Differential cross section as a function of the ratio of the two virtualities, shown as 
AQ = -logQ2/Q2 
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Figure 8: Differential cross section as a function of the acoplanarity angle between the scattered 
electrons. 
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Figure 9: Differential cross section as a function of the acolinearity $ between the scattered electrons. 
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Figure 10: Differential cross section as a function of the mass of the hadronic system W- 



77 



Oh 



T3 



25 



-§ 20 



15 



10 



— ALEPH DATA 

- - - PYTHIA 

- - PHOT02 

— NLO QCD 




%2 



Figure 1 1 : Differential cross section as a function of the deep inelastic scattering variable x% 
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Figure 12: Differential cross section as a function of the deep inelastic scattering variable yi 
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Figure 13: Differential cross section as a function of Y 
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Figure 14: The cr 77 cross section as a function of Y in comparison with LO BFKL and NLO BFKL, 
the range giving the uncertainty estimated from the variation of the Regge scale parameter from Q 2 to 
10Q 2 and from Q 2 to 4Q 2 , respectively. 
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Appendix 



Table 2: Differential cross section as a function of the relative energy of the electrons ei, 2 = E ^' 2 



e l,2 


bin 




da/de lt2 [pb] 










data 


PYTHIA 


PHOT02 


QCD 


0.30 - 


- 0.40 


2.85 


± 0.47 ± 1.03 


1.37 


1.51 


1.18 


0.40 - 


- 0.50 


2.78 


± 0.39 ± 0.30 


1.69 


1.74 


1.50 


0.50 - 


- 0.60 


3.01 


± 0.37 ± 0.32 


1.98 


2.14 


1.69 


0.60 - 


- 0.70 


3.12 


± 0.36 ± 0.10 


2.63 


2.78 


2.18 


0.70 - 


- 0.80 


4.36 


± 0.36 ± 0.21 


3.72 


3.94 


3.20 


0.80 - 


- 0.90 


5.19 


± 0.36 ± 0.27 


6.35 


6.51 


5.31 


0.90 - 


- 0.95 


11.51 


± 0.78 ± 0.40 


11.86 


12.06 


9.76 


0.95 - 


- 1.00 


27.24 ± 1.75 ± 2.99 


31.33 


29.54 


25.02 








X 2 / NDF 


23.3 / 7 


20.2 / 7 


34.6 / 8 



Table 3: Differential cross section as a function of the polar angle 9 of the scattered electrons. 



e b 


in 


da/d9 [pb/rad] 






[ra 


d] 


data 


PYTHIA 


PHOT02 


QCD 


0.035 - 


0.040 


110.11 ± 9.63 ± 3.92 


129.33 


132.93 


115.79 


0.040 - 


0.045 


122.18 ± 9.94 ± 3.01 


101.10 


104.87 


89.68 


0.045 - 


0.050 


94.71 ± 8.01 ± 2.47 


81.76 


86.08 


69.86 


0.050 - 


0.055 


81.16 ± 7.08 ± 2.73 


66.91 


69.70 


57.97 


0.060 - 


0.070 


40.57 ± 3.63 ± 1.65 


43.84 


44.06 


35.49 


0.070 - 


0.090 


28.08 ± 2.14 ± 0.92 


28.37 


28.03 


22.92 


0.090 - 


0.120 


15.81 ± 1.29 ± 0.55 


15.46 


14.38 


11.44 


0.120- 


0.155 


6.49 ± 0.84 ± 0.27 


8.30 


7.03 


5.44 






X 2 / NDF 


18.4 / 7 


13.1 / 7 


45.8 / 8 
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Table 4: Differential cross section as a function of the virtuality Q\ 2 of the photons. 



Q 2 bin 
[GeV] 






data 




da/dQf 


[pb/GeV 2 

P VTITT A 

111 Xll/l 


rnU L \JZ 




2.0 


6.0 


n m ^ 

U.UIOO 


-L 


n nnf!9 

u.uuuz 


-L 


u.uuou 


u.uuoy 


U.UIOI 


U.UUOJ 


6.0 


- 10.0 


0.0852 


± 


0.0112 


± 


0.0059 


0.0461 


0.0638 


0.0401 


10.0 


- 20.0 


0.1290 


± 


0.0070 


± 


0.0052 


0.1295 


0.1462 


0.1160 


20.0 


- 30.0 


0.0785 


± 


0.0056 


± 


0.0161 


0.0691 


0.0682 


0.0606 


30.0 


- 40.0 


0.0415 


± 


0.0081 


± 


0.0073 


0.0414 


0.0393 


0.0331 


40.0 


- 50.0 


0.0263 


± 


0.0030 


± 


0.0013 


0.0278 


0.0256 


0.0216 


50.0 


- 70.0 


0.0160 


± 


0.0016 


± 


0.0004 


0.0174 


0.0165 


0.0137 


70.0 


- 100.0 


0.0111 


± 


0.0011 


± 


0.0008 


0.0100 


0.0084 


0.0073 


100.0 


- 150.0 


0.0038 


± 


0.0005 


± 


0.0003 


0.0051 


0.0038 


0.0033 


150.0 


- 200.0 


0.0010 


± 


0.0003 


± 


0.0002 


0.0022 


0.0013 


0.0013 








X 


2 / NDF 




27.0 / 9 


12.0 / 9 


30.3 / 10 



Table 5: Differential cross section as a function of the ratio of the virtualities of the two photons 

Ql 

of- 

AQ bin da/dAQ [pb] 



AQ = -log^ 



[GeV] 




data 




PYTHIA 


PHOT02 


QCD 


-3.0 - -1.5 


0.096 


± 0.021 ± 


0.007 


0.149 


0.120 


0.096 


-1.5 - -1.0 


0.394 


± 0.057 ± 


0.016 


0.374 


0.357 


0.292 


-1.0 - -0.5 


0.631 


± 0.076 ± 


0.025 


0.519 


0.546 


0.464 


-0.5 - -0.2 


0.686 


± 0.098 ± 


0.042 


0.619 


0.661 


0.565 


-0.2 - 0.0 


0.749 


± 0.122 ± 


0.034 


0.623 


0.755 


0.611 


0.0 - 0.2 


0.664 


± 0.127 ± 


0.044 


0.638 


0.763 


0.633 


0.2 - 0.5 


0.736 


± 0.105 ± 


0.043 


0.612 


0.709 


0.566 


0.5 - 1.0 


0.620 


± 0.083 ± 


0.035 


0.524 


0.538 


0.451 


1.0 - 1.5 


0.403 


± 0.062 ± 


0.023 


0.383 


0.354 


0.276 


1.5 - 3.0 


0.122 


± 0.021 ± 


0.005 


0.147 


0.120 


0.095 






X 2 / NDF 




13.1 / 9 


4.7 / 9 


20.9 / 10 
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Table 6: Differential cross section as a function of the acolinearity $ of the scattered electron and the 
scattered positron. 



$ bin 








da/d® 


[pb/deg 


ree] 




[de 


gree] 






data 




PYTHIA 


PHOT02 


QCD 


166.0 


- 170.0 


0.014 


± 


0.004 ± 


0.001 


0.017 


0.017 


0.015 


170.0 


- 173.0 


0.088 


± 


0.010 ± 


0.002 


0.090 


0.102 


0.081 


173.0 


- 174.0 


0.189 


± 


0.029 ± 


0.013 


0.184 


0.230 


0.172 


174.0 


- 175.0 


0.307 


± 


0.032 ± 


0.021 


0.257 


0.316 


0.251 


175.0 


- 176.0 


0.295 


± 


0.033 ± 


0.025 


0.301 


0.376 


0.281 


176.0 


- 177.0 


0.308 


± 


0.034 ± 


0.034 


0.271 


0.293 


0.223 


177.0 


- 178.0 


0.220 


± 


0.030 ± 


0.014 


0.266 


0.217 


0.184 


178.0 


- 179.0 


0.200 


± 


0.031 ± 


0.016 


0.229 


0.120 


0.138 


179.0 


- 179.5 


0.192 


± 


0.059 ± 


0.070 


0.170 


0.066 


0.093 








X 2 / NDF 




5.9 / 8 


15.6 / 8 


12.1 / 9 



Table 7: Differential cross section as a function of the acoplanarity angle A$ of the scattered electron 
and the scattered positron. 



Acf) bin 






da/dA(f> [pb/degree] 




[degree] 




data 




PYTHIA 


PHOT02 


QCD 


0.00 - 0.31 


0.45 


± 0.07 ± 


0.02 


0.42 


0.27 


0.43 


0.31 - 0.63 


0.49 


± 0.07 ± 


0.02 


0.44 


0.36 


0.45 


0.63 - 0.94 


0.58 


± 0.08 ± 


0.06 


0.45 


0.77 


0.48 


0.94 - 1.26 


0.54 


± 0.08 ± 


0.07 


0.49 


0.85 


0.55 


1.26 - 1.57 


0.53 


± 0.08 ± 


0.06 


0.53 


0.79 


0.60 


1.57 - 1.88 


0.62 


± 0.08 ± 


0.05 


0.61 


0.82 


0.58 


1.88 - 2.20 


0.82 


± 0.09 ± 


0.03 


0.68 


0.88 


0.56 


2.20 - 2.51 


0.58 


± 0.09 ± 


0.03 


0.76 


0.83 


0.53 


2.51 - 2.83 


0.83 


± 0.10 ± 


0.08 


0.89 


0.41 


0.50 


2.83 - vr 


0.88 


± 0.13 ± 


0.14 


0.99 


0.29 


0.48 






X 2 / NDF 




8.7 / 9 


60.7 / 9 


19.5 / 10 
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Table 8: Differential cross section as a function of the mass W 77 of the hadronic system. 



W 77 bin 


da/dWj-f 


[pb/GeV] 






[GeV] 


data 


PYTHIA 


PHOT02 


QCD 


3.0 - 5.0 


0.1409 ± 0.0164 ± 0.0257 


0.1499 


0.1307 


0.0996 


5.0 - 10.0 


0.1263 ± 0.0095 ± 0.0033 


0.1334 


0.1380 


0.1019 


10.0 - 15.0 


0.0763 ± 0.0072 ± 0.0032 


0.0787 


0.0814 


0.0665 


15.0 - 35.0 


0.0276 ± 0.0023 ± 0.0023 


0.0246 


0.0244 


0.0205 


35.0 - 50.0 


0.0102 ± 0.0020 ± 0.0021 


0.0050 


0.0052 


0.0043 


50.0 - 80.0 


0.0005 tomol ± °- 0020 


0.0012 


0.0013 


0.0011 



X 2 / NDF 5.3 / 5 6.4/5 18.0 / 6 



Table 9: Differential cross section as a function of the Bjorken variable xi t 2- 





bin 


data 


da/dxi [pb] 
PYTHIA 


PHOT02 


QCD 


0.00 


- 0.05 


17.71 ± 1.38 ± 


2.00 


14.11 


15.88 


12.17 


0.05 


- 0.10 


14.48 ± 1.05 ± 


1.13 


13.51 


13.70 


10.84 


0.10 


- 0.15 


11.97 ± 0.95 ± 


0.77 


11.12 


11.00 


8.99 


0.15 


- 0.20 


8.20 ± 0.79 ± 


0.44 


8.67 


8.62 


7.02 


0.20 


- 0.25 


6.62 ± 0.72 ± 


0.56 


7.03 


6.98 


5.61 


0.25 


- 0.30 


7.31 ± 0.70 ± 


0.43 


5.40 


5.50 


4.66 


0.30 


- 0.40 


3.71 ± 0.36 ± 


0.22 


3.68 


3.96 


3.45 


0.40 


- 0.50 


2.19 ± 0.28 ± 


0.14 


2.25 


2.31 


2.07 


0.50 


- 1.00 


0.42 ± 0.06 ± 


0.09 


0.70 


0.46 


0.46 






X 2 / NDF 




16.0 / 8 


7.3 / 8 


30.6 / 9 



Table 10: Differential cross section as a function of the Bjorken variable 2/1,2- 



bin 



da/dyi [pb] 





data 




PYTHIA 


PHOT02 


QCD 


0.000 


- 0.025 


51.28 ± 5.54 ± 


7.34 


40.37 


37.49 


31.79 


0.025 


- 0.050 


19.41 ± 1.99 ± 


2.07 


21.48 


20.73 


17.61 


0.050 


- 0.100 


11.68 ± 0.79 ± 


0.49 


12.02 


12.27 


9.93 


0.100 


- 0.150 


5.16 ± 0.50 ± 


0.35 


7.54 


7.55 


6.19 


0.150 


- 0.200 


5.39 ± 0.52 ± 


0.21 


5.31 


5.61 


4.49 


0.200 


- 0.300 


4.43 ± 0.36 ± 


0.22 


3.74 


3.94 


3.21 


0.300 


- 0.400 


3.15 ± 0.36 ± 


0.09 


2.64 


2.79 


2.20 


0.400 


- 0.600 


2.85 ± 0.27 ± 


0.26 


1.84 


1.95 


1.60 


0.600 


- 0.800 


1.41 ± 0.23 ± 


0.50 


0.69 


0.76 


0.59 






X 2 / NDF 




31.1 / 8 


28.1 / 8 


42.4 / 9 
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Table 11: Differential cross section as a function of Y. 



Y bin 



da/dY [pb] 



data 



QCD 



-1.0 


- 0.0 


0.224 ± 0.028 ± 0.031 


0.270 


0.223 


0.195 


0.0 


- 1.0 


0.437 ± 0.039 ± 0.015 


0.483 


0.493 


0.381 


1.0 


- 2.0 


0.542 ± 0.045 ± 0.022 


0.523 


0.545 


0.425 


2.0 


- 3.0 


0.422 ± 0.040 ± 0.030 


0.352 


0.370 


0.312 


3.0 


- 4.0 


0.227 ± 0.035 ± 0.017 


0.171 


0.194 


0.149 


4.0 


- 5.0 


0.117 ± 0.026 ± 0.024 


0.062 


0.079 


0.065 


5.0 


- 6.0 


0.026 tomo ± 0.013 


0.019 


0.025 


0.012 


6.0 


- 7.0 


0.008 t o ol 8 ± 0.005 


0.004 


0.009 


0.004 



x7 NDF 



9.3 / 7 5.2/7 19.3 / 
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